Renormalization Group Theory for the Imbalanced Fermi Gas 
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We formulate a wilsonian renormalization group theory for the imbalanced Fermi gas. The theory 
is able to recover quantitatively well-established results in both the weak-coupling and the strong- 
coupling (unitarity) limit. We determine for the latter case the line of second-order phase transitions 
of the imbalanced Fermi gas and in particular the location of the tricritical point. We obtain good 
agreement with the recent experiments of Y. Shin et al. [Nature 451, 689 (2008)]. 
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Introduction. The amazing experimental con- 

trol in the manipulation of degenerate Fermi mixtures, 
has led for the balanced mixture to an accurate study 
of the crossover between a Bardeen-Cooper-Schrieffer 
(BCS) superfluid and a Bose-Einstein condensate (BEC) 
of diatomic molecules. Particularly interesting is the 
strongly-interacting regime, where the scattering length 
of the interaction becomes much larger than the average 
interatomic distance. In this so-called unitarity limit, 
experiments have revealed that the superfluid state is re- 
markably stable and has a record-high critical tempera- 
ture of about one tenth of the Fermi energy [l[ . Theoreti- 
cally, the unitarity limit is extremely challenging, because 
there is no rigorous basis for perturbation theory due to 
the lack of a small parameter. As a result, mean-field the- 
ory is only useful for understanding the relevant physics 
qualitatively, but cannot be trusted quantitatively. In or- 
der to get accurate results, more sophisticated theoretical 
methods have to be invoked. 

An important example is using quantum Monte-Carlo 
techniques, which can provide exact results about the 
strongly- interacting regime @, H, 0, HI' but offer less 
physical insight than analytic methods. Therefore, sev- 
eral other approaches have been developed to improve 
on mean-field theory. Examples are theories incorporat- 
ing Gaussian fluctuations [q, 0, H M, E3], e expansion 
11 1, l/N expansion fl2j fl3l , and the functional renor- 
malization group (RG) [hj, |l5j. In this Letter, we for- 
mulate a so-called wilsonian RG to study the strongly- 
interacting atomic Fermi mixture with a population im- 
balance. The intuitively appealing wilsonian approach, 
which has been extremely successful in the study of crit- 
ical phenomena (l6| . is based on systematically integrat- 
ing out short-wavelength degrees of freedom, which then 
renormalize the coupling constants in the effective action 
for the long-wavelength degrees of freedom. For fermions, 
the excitations of lowest energy lie near the Fermi level, 
which is therefore the natural end point for a renormal- 
ization group flow 17 1 - A notorious problem for inter- 



We show, however, how to obtain RG equations that au- 
tomatically flow to the final value of the renormalized 
Fermi level. 

The unitary, two-component Fermi mixture with an 
unequal number of particles in each spin state is a topic 
of great interest in atomic physics, condensed matter, 
nuclear matter, and astroparticle physics. The landmark 
atomic-physics experiments exploring this system, per- 
formed at MIT by Zwicrlcin et al. [2fJ and at Rice Uni- 



versity by Partridge et al. 2l|, induced a large amount 



acting fermions is that under renormalization the Fermi 
level also flows to an a priori unknown value, making 
the wilsonian RG difficult to perform in practice 18, 1||. 



of activity, caused by an intriguing mix of mutual consis- 
tent and contradictory results. In summary, both exper- 
iments observed no oscillating order parameter, so that 
the Fulde-Ferrell and Larkin-Ovchinnikov phases do not 
seem to play a role in the unitarity limit. Therefore, the 
experiments are consistent with a phase diagram includ- 
ing both second-order and first-order phase transitions 
between the superfluid (BCS or Sarma) phase and the 
normal phase, that are connected by a tricritical point 
[ij, l22j]. However, as a function of population imbalance 
Zwierlein et al. obtain a critical imbalance at which the 
trapped Fermi gas becomes fully normal, whereas Par- 
tridge et al. observe a superfluid core up to their highest 
imbalances. Although this contradictory result is still not 
completely understood, more recent work implies that 
the data of Zwierlein et al. is consistent with the local- 
density approximation, whereas the experiments of Par- 
tridge et al. explore physics beyond this approximation, 
possibly due to the smaller number of particles and the 
more extreme aspect ratio of the trap 2^, 24 1 . 

Since the validity of the local-density approximation 
implies that the Fermi mixture can be seen as being lo- 
cally homogeneous, the MIT group is in the unique posi- 
tion to experimentally map out the homogeneous phase 
diagram by performing local measurements in the trap. 
Most recently, this important experiment was performed 
by Shin et al. [251 ] . obtaining for the homogeneous tri- 
critical point in the unitarity limit P C 3 = 0.20(5) and 
T C 3 = 0.07(2) Tpj, with P the local polarization given 
by P = (nf — «x)/( n T ~b n l)i n o t ne density of atoms in 
spin state \a), T the temperature, and ep a = k^,T-p a = 
(67r 2 n CT ) 2 / 3 7i 2 /2m the Fermi energies with m the atomic 
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mass. So far, there has not been an accurate calculation 
for this homogeneous tricritical point. In this Letter, we 
determine it to lie at P C 3 = 0.24 and T c3 = 0.06 Tpj, in 
good agreement with the experiment by Shin et al. 

Wilsonian renormalization. - The central idea of 
wilsonian renormalization is to subsequently integrate 
out degrees of freedom in shells at high momenta A of 
infinitesimal width dA and absorb the result of the inte- 
grations into various coupling constants, which are there- 
fore said to flow. First, we calculate the Feynman dia- 
grams renormalizing the coupling constants of interest, 
while keeping the integration over the internal momenta 
restricted to the considered high-momentum shell. Only 
one-loop diagrams contribute to the flow, because the 
thickness of the momentum shell is infinitesimal and each 
loop introduces a factor dA. In order to obtain the ex- 
act partition sum, it is then needed to consider an infi- 
nite number of coupling constants. Although this is not 
possible in practice, the RG is still able to distinguish 
between the relevance of the various coupling constants, 
such that a carefully selected set of them already leads 
to highly accurate results. 

Consider the action of an interacting Fermi mixture 

SW,<j>] = ^2 Ck,nH^«+ e k-/v)^k,n (1) 
k,n,o" 

+ *- r q , m 0* q _ k / m _„/0J k / „/(/)J,q-k,m-n0Tk,n j 
k,k',q 

n,n ,m 

with u) n the odd fermionic Matsubara frequencies, e k = 
h 2 k 2 /2m the kinetic energy, \i a the chemical potentials, 
(3 = l/k B T, V the volume, r q m the interaction ver- 
tex and 4> a k,n the fermionic fields corresponding to an- 
nihilation of a particle with spin er, momentum k and 
frequency ui n . In Fig. 1, we have drawn the Feynman 
diagrams renormalizing fi a and r q m . To start with a 
simple wilsonian RG, we take the interaction vertex to 
be frequency and momentum independent. If we then 
consider only the three coupling constants \i a and To^, 
we find 



dA 
~dA 



2ir 2 



1 — Ni- — N\ 



Ni 



2(e A - m) 



2/; 



A 2 

2^ IV 



(2) 
(3) 



with [i = [hi + /ijJ/2, h — (/x-f — ^|)/2 and the Fermi 
distribution N a = l/{exp[/3(eA — )] + !}• These expres- 
sions are readily obtained from the diagrams in Fig. 1 
by setting all external frequencies and momenta equal to 
zero and by performing in each loop the full Matsubara 
sum over internal frequencies, while integrating the in- 
ternal momenta over the infinitesimal shell dA. The first 
term in Eq. ^ corresponds to the ladder diagram and 
describes the scattering between particles. The second 
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FIG. 1: Feynman diagrams renormalizing a) the chemical po- 
tentials and b) the interatomic interaction. 



term corresponds to the bubble diagram and describes 
screening of the interaction by particle- hole excitations. 
Also note that due to the coupling of the differential equa- 
tions for /v and Tq g, we automatically generate an infi- 
nite number of Feynman diagrams, showing the nonper- 
turbative nature of the RG. 

However, when the Fermi mixture is critical, the in- 
verse many-body vertex Tq q flows to zero according to 
the Thouless criterion and the chemical potentials in Eq. 
([3]) diverge, which is unphysical. To go beyond this sim- 
ple RG and calculate critical properties realistically, we 
need to take the frequency and momentum dependence 
of the interaction vertex into account, which are gener- 
ated by the ladder and the bubble diagrams. The lad- 
der diagram depends only on the external center-of-mass 
coordinates q and uj m , and its contribution to the renor- 
malization of is given by 



E(q 2 ,iu m ) 



r/q' 



, l-iV T (e q ,)-iV|(e q - q 



-ihcdr, 



(4) 



where during integration both q' and q — q' have to re- 
main in the infinitesimal shell dA. Since the ladder dia- 
gram is already present at the two-body level, it is most 
important for scattering properties. Therefore, the inter- 
action vertex is mainly dependent on the center-of-mass 
coordinates and we neglect the dependence of the vertex 
on other frequencies and momenta. The way to treat the 
external frequency and momentum dependence in wilso- 
nian RG, is by expanding the (inverse) interaction in the 



following way: 



r -i 



= r, 



O,0 



-l q 2 



Z,,,itiu) m - The 



flow equations for the additional coupling constants Z q 1 
and Z~ x are then obtained by considering the derivatives 
d q 2~,(q 2 ,u)\ q=UJ= o and d ul E(q 2 ,uj)\ q=UJ=0 . 

Extreme imbalance. — First, we apply the RG to one 
spin-down particle in a Fermi sea of spin-up particles at 
zero temperature in the unitarity limit. The full equa- 
tion of state for the normal state of a strongly-interacting 
Fermi mixture was obtained at zero temperature using 
Monte-Carlo techniques 0]. The most important feature 
of this equation is a so-called mean-field shift, caused by 
the strong interactions and characterized by a parameter 
A, which describes the self-energy of a single spin-down 
particle in a sea of spin- up particles 0, [26J. In this case, 
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the RG equations are simplified, because N± can be set 
to zero and thus \i\ is not renormalized. Next, we have to 
incorporate the momentum and frequency dependence of 
the interaction in the one-loop Feynman diagram for the 
renormalization of Hi- In this particular case, the exter- 
nal frequency dependence of the ladder diagram can be 
taken into account exactly, since the one-loop Matsub- 
ara sum simply leads to the substitution %hio m — ► e q — H\ 
in Eq. ((4]) [26(. The external momentum dependence is 
accounted for by the coupling Z~ x , giving 



(5) 
(6) 

(7) 
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where we note that these equations only have poles for 
positive values of fit. Since this will not occur, we can 
simply use A(Z) = A e - ' to integrate out all momen- 
tum shells [27|. We then obtain a system of three cou- 
pled ordinary differential equations in I, which are very 
easily solved numerically. If we take as an initial con- 
dition r o(0) = -m(it + 2\a\A ) / An 2 \a\h 2 for a nega- 
tive scattering length a, we automatically incorporate 
the relevant two-body physics exactly into our theory 
and also eliminate all dependence on the high-momentum 
cut-off Ao 27]. The unitarity limit is then given by 
Tq q(0) = — mAo/2Tr 2 h 2 . The other initial conditions are 
AH(0) = /ij and 2'~ 1 (0) = 0, since the interaction starts 
out as being momentum independent. Note that in this 
calculation /Lt^(0) = Hi is indeed negative and increases 
during the flow due to the strong attractive interactions. 
The quantum phase transition from a zero density to a 
nonzero density of spin-down particles occurs for the ini- 
tial value hi that at the end of the flow precisely leads to 
Hi(oo) = 0. This happens when hi — — 0.598/Zf, yielding 
A = 0.997 in very good agreement with the Monte Carlo 
result A = 0.97(2) ||. This calculation also shows that 
it is crucial to let the chemical potential flow. 

Phase diagram. — Next, we turn to our main topic, 
namely the critical properties of the strongly-interacting 
Fermi mixture and the calculation of the tricritical point 
in the phase diagram. Since it is not exact to make the 
substitution Hlu — > e q — h-<j a t nonzero temperatures, 
we take the frequency dependence of the ladder diagram 
into account through the renormalization of the coupling 
Zj 1 . While the flow of Y~^\ is still given by Eq. ([2]), the 
expressions for the flow of h<j and Z~ x become 
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dA 



A 2 



N-„ + N B 



2ir 2 
A 2 1 - N- 



1 0,0 



Z^A 2 



Zu (e A - H-a) 



2tt 2 A{e K -H? 



(8) 
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with A^b = l/{exp[/3^(-r oo + Z- X A 2 )} - 1} coming 
from the bosonic frequency dependence of the interac- 
tion. A more cumbersome expression holds for Z~ x . The 
initial conditions are the same as for the extremely imbal- 
anced case with in addition /Ltf(0) = Hi an d Zj 1 (0) = 0. 
As mentioned before, the critical condition is that the 
fully renormalized vertex Tq (oo), which can be seen as 
the inverse many-body T-matrix at zero external momen- 
tum and frequency, goes to zero. Physically, this implies 
that a many-body bound-state is entering the system. 
From Eq. ©, we see that incorporating the coupling 
constants Z~ x and Z~ x , thereby taking the dependence 
of the interaction on the center-of-mass momentum and 
frequency into account, is crucial to solve the previously 
mentioned problem of the diverging chemical potential. 

The only pole left in our set of RG equations is the 
average Fermi level h — if 1 ] + which is there- 

fore the natural end point of our RG. However, this 
Fermi level is shifting due to the renormalization of the 
individual chemical potentials. This problem is conve- 
niently solved by integrating out all momentum shells 
with the following procedure. First, we start at a high 
momentum cutoff Ao and flow to a momentum A' at 
roughly two times the average Fermi momentum, with 
k-pa = \/2meF (7 /h~. This integrates out the high-energy 
two-body physics, but hardly affects the chemical po- 
tentials. Then, we start integrating out the rest of 
the momentum shells symmetrically with respect to the 
flowing average Fermi level. This is achieved by using 
hA + (l)/V2m = (hA' /^/2m - ^fi)e~ l + y^MO and by 
HA-(l)/y / 2m = —y/Jle~ l + \J Note that as desired 
A+(Z) starts at A' Q and automatically flows from above to 
y/2mH{oo)/h, whereas A_(Z) starts at and automati- 
cally flows from below to ^/2m^(oo)/?i. 

We first apply the above procedure to study the equal 
density case, i.e., h — > 0, as a function of negative scatter- 
ing length a. The scattering length enters the calculation 
through the initial condition of Tq . To express our re- 
sults in terms of the Fermi energy ep = cfo-, we calculate 
the densities of atoms with the flow equation dn a /dA = 
A 2 N rT /2ir 2 . In the weak-coupling limit, a — > 0~, the 
chemical potentials hardly renormalize, so that only Eq. 
(|2"|) is relevant. The critical temperature becomes expo- 
nentially small, which allows us to solve Eq. © exactly 
with the result k^T c = 8epe 7-3 exp{— 7r/2fcF|a|}/7r and 7 
Euler's constant. Compared to the standard BCS-result 
we have an extra factor of 1/e, coming from the screening 
effect of the bubble diagram that is not present in BCS 
theory. It is to be compared with the so-called Gor'kov 
correction, that reduces the critical temperature by a fac- 
tor of 2.2 in the weak-coupling BCS-limit [II]. The dif- 
ference is due to the fact that we have only allowed for a 
nonzero center-of-mass momentum. This approximation 
is actually most appropriate in the unitarity limit and 
expected to be less accurate for weak coupling. 
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FIG. 2: (Color online) The phase diagram of the homoge- 
neous two-component Fermi mixture in the unitarity limit, 
consisting of the superfluid Sarma (S) and BCS phases, the 
normal phase (N) and the forbidden region (FR). The solid 
black line is the result of our RG calculations. The Monte- 
Carlo result of Lobo et al. [fj, which is recovered by our RG, 
is indicated by a cross. The open circles and squares are data 
along the phase boundaries from the experiment of Shin et al. 
[25l | . The dashed lines are only guides to the eye. Also shown 
is the Feynman diagram determining the tricritical point. 



At larger values of |a|, the flow of the chemical po- 
tential becomes important and we obtain higher critical 
temperatures. In the unitarity limit, when a diverges, we 
obtain T c = 0.13Tp and ^(T c ) = 0.55cf in good agree- 
ment with the Monte-Carlo results T c = 0.152(7)T F and 
M (T C ) = 0.493(14)e F [4]. With our RG approach, we are 
in the unique position to accurately calculate the critical 
temperature as a function of polarization P and compare 
with the recent experiment of Shin et al. The result is 
shown in Fig. 2. The inset of this figure shows the one- 
loop diagram determining the position of the tricritical 
point. If it changes sign, then the fourth-order coeffi- 
cient in the Landau theory for the superfluid phase tran- 
sition changes sign and the nature of the phase transi- 
tion changes from second order to first order. This yields 
P C 3 = 0.24 and T c $ — 0.06 Tpj in good agreement with 
the experimental data. Our previous confirmation of the 
Monte-Carlo equation of state at T — implies that we 
also agree with the prediction of a quantum phase transi- 
tion from the superfluid to the normal phase at a critical 
imbalance of P c = 0.39 [5]. Note that up to now, all 
theoretical predictions for the location of the tricritical 
point do not fit on the scale of Fig. 2. However, our 
calculations find good agreement with the experiments 
of Shin et al. in all limits. 

Near the second-order phase boundary, the BCS or- 
der parameter |A| becomes arbitrarily small. Since at 
nonzero polarization we have that /i(oo) > 0, it immedi- 



ately follows that h(oo) > |A|. This means that the nor- 
mal gas is unstable towards the so-called Sarma phase, 
which is a polarized superfluid with a ga ples s excitation 
spectrum for the majority spin-species [221 ]. However, 
the present RG is not suitable to calculate the full extent 
of the Sarma phase in the phase diagram or the precise 
shape of the forbidden region, because this requires a RG 
for the superfluid phase. The corresponding calculations 
are more involved than the present RG for the normal 
phase and is work in progress. 
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